Exotic quantum phase transitions in a Bose-Einstein condensate coupled to an optical cavity 



00 

o 

O 



-f— > 



(N 
> 

t> 

On 

m 

O 

oo 
o 



X 



Gang Chen, 1 - 2 Xiaoguang Wang, 3 J. -Q. Liang, 2 and Z. D. Wang ''El 

'Department of Physics and Center of Theoretical and Computational Physics, 
The University of Hong Kong, Pokfulam Road, Hong Kong, China 
2 Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China 
3 Zhejiang Institute of Modern Physics, Department of Physics, Zhejiang University, Hangzhou 310027, China 

A new extended Dicke model, which includes atom-atom interactions and a driving classical laser field, is 
established for a Bose-Einstein condensate inside an ultrahigh-finesse optical cavity. A feasible experimental 
setup with a strong atom-field coupling is proposed, where most parameters are easily controllable and thus 
the predicted second-order superradiant-normal phase transition may be detected by measuring the ground-state 
atomic population. More intriguingly, a novel second-order phase transition from the superradiant phase to the 
"Mott" phase is also revealed. In addition, a rich and exotic phase diagram is presented. 
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As is known, a trapped Bose-Einstein condensate (BEC) 
may be used to generate a macroscopic quantum object con- 
sisting of many atoms that are in the same quantum state with 
a longer lifetime and can be excited by either deforming the 
trap or varying the interactions among atoms. Thus the BEC, 
as a distinct macroscopic quantum system, plays an important 
role in the in-depth exploration of both fundamental physics 
and quantum device applications of many -body systems J2]. 
In particular, an intriguing idea to combine the cavity quan- 
tum electrodynamics (QED) with the BEC has recently at- 
tracted significant interests both theoretically and experimen- 
tally as many exotic quantum phenomena closely related to 
both light and matter at ultimate quantum levels may emerge 
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Very recently, a so-called strong coupling of a BEC to the 
quantized field of an ultrahigh-finesse optical cavity was re- 
alized experimentally |fl4|, which not only implies that a new 
challenging regime of cavity QED has been reached, where all 
atoms occupying a single mode of a matter-wave field that can 
couple identically to the photon induced by the cavity mode, 
but also opens a wider door to explore a variety of new quan- 
tum phenomena associated with the cavity-mediated many- 
body physics of quantum gas. Regrettably, the authors of 
Ref. H 1 4fl ignored the important nonlinear interactions among 
the untracold atoms that are controllable via the Feshbach 
resonance technique, while these interactions are believed to 
have also a considerable impact on physical properties of the 
BEC, leading to some exotic quantum phenomena 115 1. 

In this Letter, we establish an extended Dicke model with 
the atom-atom interactions and a driving classical laser field 
under the two-mode approximation. A feasible experimen- 
tal setup with controllable parameters including a collective 
strong atom-field coupling is proposed. We illustrate how to 
drive a well-known second-order superradiant-normal phase 
transition and how to detect it experimentally. Remarkably, 
this superradiant phase transition was predicted in quantum 
optics ma ny y ears ago, but has never been observed in ex- 



periments J3, [H \m El US EH Q. More intriguingly 



phase diagram covering phenomena from quantum optics to 
the BEC, which is attributed to the competition between the 
atom-atom and the atom-field interactions. 
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novel second-order superradiant to "Mott" phase transition is 
also revealed. In addition, we also obtain a rich and exotic 



FIG. 1: (Color online) Schematic diagram of an experimental setup 
for a BEC of 87 Rb atoms coupled to a QED cavity. The BEC with 
two levels |1) and |2) is prepared in a time-averaged, orbiting poten- 
tial magnetic trap. After moving the BEC into an ultrahigh-finesse 
optical cavity, an external controllable classical laser is applied to 
produce various transitions of the atoms between |1) and 2) states. 



Our proposed experimental setup is depicted in Fig. 1. For 
an optical cavity with length 176 /jm and the mode waist 
radius 27 /im, we may choose the parameters of the cav- 
ity (g , k, i) = 2tt x (10.6,1.3,3) MHz 0], where g is 
the maximum single atom-field coupling strength, k and 7 
are the amplitude decay rates of the excited state and the in- 
tracacity field, respectively. Such a choice implies that the 
system is in the strong coupling regime, and thus the long- 
range coherence could be well established and the quantum 
dissipation effect may be safely neglected. Based on a pair 
of coupled Gross-Pitaevskii equations for the BEC with two 
levels \F = l,m f = -1) (|1)) and \F = 2,771/ = 1) (|2)) of 
5 2 <Si/2 I23h and under the two-mode approximation, the total 



Hamiltonian for the elastic two-body collisions with the inter- 
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action potential of (5-functional type may be written as 

H = H p h + Hat-ph + Hat + Hat-cl + H a t-at 



(1) 



with H p h = ua^a (h = 1 hereafter), H at = ojic\ci + 
(uj 2 + oji2)c\c 2 , H a t-at = r/ic{cJciCi/2 + r]2clclc 2 c 2 /2 + 
XC1C1C2C2, H at -d = ilfcjCi exp(— ivut) + c\c2 exp(imt)]/2, 
and H a t-ph = A(c|c2 + c^ci)^ + a), where a is the 
annihilation operator of the cavity mode with frequency 
u>; c\ and c 2 are the annihilation boson operators for |1) 
and |2), respectively; ui = J d 3 r{<p*(r)[— V 2 /2m R + 
V(r)]<f>i(r) (I = l,2)with V(r) being a single mag- 
netic trapped potential of frequencies u>i(i = x,y,z) 
and m R being the atomic mass; u>i 2 is the atomic res- 
onance frequency; 77; = (4npi/mn) J d 3 r |0/(r)| 4 and 
X = {^pi. 2 /m R ) J d 3 r |0i (r)| 2 |0 2 (r)| 2 with pi andpi, 2 (= 
P2,i) being the intraspecies and the interspecies s— wave scat- 
tering lengths, respectively; = 2Slo / d 3 r<fi 2 '(r)(f>i (r) with 
Oo being the Rabi frequency for the introduced classical laser 
with a driving frequency m; and A = g J d 3 r<f>2(r)c/)i(r) = 
g J d 3 r(fil(r)(f)2(r) with g being a interaction constant be- 
tween the atom and the photon II2411 . 

Under a unitary transformation U = cxp(—iwJ z t) with 
the condition w « u> and using the Schwinger relations 
J x = (4 c i + c I c 2)/2, J y = (c[c 2 - c\ci)/2i, and J z = 
(c\ci - c\c 2 )/2 with the Casimir invariant J 2 = N{N/2 + 
l)/2, Hamiltonian (1) can approximately be rewritten as 



y/N 



u J z + QJ x + —J 2 (2) 



H = Lua^a 



in the rotating frame, where A = 2Ay TV denotes a collective 
coupling strength, v = N[(r)i + f]2)l2 — x) describes the 
atom-atom interactions including the repulsive (v > 0) and 
attractive (v < 0) interactions, and ujq = u>2 — ^i + (N — 
l)( ? /2 — ?7i)/2 + A with A = W12 — txj being the detuning. 
For a single trapped potential, we have 0J2 = u>i and consider 
only the case of p\ = p 2 , which has the advantages that it 
reduces the effects of fluctuations in the total atomic number 
and ensures a large spatial overlap of different components 
of the condensate wavefunction. Thus, the parameters v and 
u>o can further be reduced to v = N(r/i — x) an d = A. 
Eq. (2) is a key result, which describes the collective dynamics 
for the composite system and has a rich phase diagram. Here 
we refer this equation to as an extended Dicke model since 
it contains the extra laser field term (the 4-th one) and atom- 
atom interaction term (the 5-th one) in comparison with the 
standard Dicke model and its generalized version iflH Eoll . 

A distinct property of Hamiltonian (2) lies in that all pa- 
rameters can be controlled independently. For example, the 
effective coupling strength A can be manipulated by a stan- 
dard technique. The effective Rabi frequency Q, and the de- 
tuning A depend on the experimentally controllable classical 
laser, and especially, the detuning A can vary continuously 
from the red (A < 0) to the blue (A > 0) detunings. The 
parameter v ranging from the positive to the negative is deter- 



mined by the s— wave scattering lengths via Feshbach reso- 
nance technique |15j] . For v = (pi.2 = Pi) and fi = w = 0, 
Hamiltonian (2) is reduced to a standard Dicke model with a 
second-order superradiant phase transition at the critical point 



)na-oraer superradiant pnase tr 
A c = ^/UJuJE G3 E3 H H El IHl |2J . It should be no- 
ticed that this important prediction has never been observed 
in experiments. The main difficulties are likely (i) all atoms 
can hardly interact identically with the same quantum field; 
(ii) the frequencies u> and loq typically exceed the coupling 
strength A by many orders of magnitude; (iii) it is hard to con- 
trol the parameters as demanded. However, in our proposal, 
these difficulties could be completely overcome by using the 
currently available experimental techniques of BEC, as will 
readily be seen below. 

To explore quantum phases and their transitions, we 
now investigate the ground-state properties of Hamiltonian 
(2), which can approximately be dealt with by using the 
Holstein-Primakoff transformation, Jj = b' ^ N — b^b, J_ = 
ViV - Wbb and J z = {tfb - N/2) with [b, 6 f ] = 1. Here 
we introduce two shifting boson operators c) = a) + \/~N a 
and = b^ — yN {3 with auxiliary parameters a and /3 to 
describe the collective behaviors of both the atoms and the 
photon H IH Q. With the help of the boson ex- 
pansion method, the scaled ground-state energy is given by 

„2 



E (a,h)/N = 
n(h 2 - 1/2) + 
(1/2 < h 2 < 
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2Xa{h 2 - 1/2) 
- h 2 ) with h\ft~- 
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h 2 = [3 2 - 1/2 
The critical points can be determined 
from the equilibrium condition d[Eo(a, h)/N]/da = and 
d[Eo(a, h) / N] / dhx dh/ df3=Q, which leads to two equations: 
a = A(?7 2 - l)/2w(?7 2 + 1) and 



2(u + v)r)(l - rf) + 2fi»/(l + rf) + A(l - rf) = 0, (3) 

where rj — h/yl — h 2 and u = \ 2 /uj are introduced as 
new parameters for convenience. The coefficient (u + v) de- 
scribes the intrinsic competition between the atom-atom and 
the atom-field interactions and gives rise to some exotic phase 
transitions predicted in the following. 

Equation (3) contains the basic information of quantum 
phases and their transition. As a benchmark, we first address 
the simplest case that there is no nonlinear interaction among 
atoms, namely, v = (p± = p\.2)- Fig. 2 shows the scaled 
ground-state energy Eq/N and atomic population (or equiva- 
lently "magnetization") m/N as a function of the detuning A 
for different Rabi frequencies (£1). It can be seen clearly that 
in the limit f2 — > 0, this system exhibits collective excitations 
of both the atom and the field with macroscopic occupations 
(i.e., \m/N\ < 1 and (a'a) > 0) for —it < A < u, whereas 
there are no such excitations for A > u and A < — u (the 
solid black line). This interesting behavior typically shows the 
second-order superradiant phase transition in quantum o ptic s 
with the critical point A c = ±u lfl6l \vj [HI [l9i |20j |2lj Wa . 
Moreover, here we may achieve the condition that the order 
of magnitude of A is the same as that of \/ujA by controlling 
the detuning of the classical laser. By controlling Sl/s C 1 
and evaluating a partial derivative of m with respect to A ( 
or w), if a peak is detected in the derivative, which becomes 
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FIG. 2: (Color online) The scaled ground-state energy Eo/N and 
atomic population m/N (Insert part) versus the detuning A for dif- 
ferent Rabi frequencies (Q.) at v = 0. 



sharper and shaper if CI becomes smaller and smaller, a sec- 
ond order superradiant phase transition at CI = is signa- 
tured, even though the transition disappears at a finite Rabi 
frequency CI. In view of this, our proposed composite system 
with the controllable classical laser is a promising candidate 
for exploring cavity-induced superradiant phase transition by 
measuring the ground-state atomic population via the resonant 
absorption imaging 1 24] . 
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FIG. 3: (Color online) A zero- temperature phase diagram of the de- 
tuning A and the atom-atom interaction strength v with a Rabi fre- 
quency CI. The blue line is determined by A = ±[(u + v) 2 ^ 3 — 
fi 2 ' 3 ] 3 ' 2 . Two interesting second-order phase transitions occur 
when crossing the red and blue lines. Note that this phase dia- 
gram is symmetric with respect to A < 0, simply because the 
Hamiltonian (2) is invariant under the transformation A = —A 
and J z = —Jz- The lower Table denotes three different quan- 
tum phases. Here, "SF"="Superfluid" phase, "M"="Mott" phase, 
SR=Superradiant phase, P=phase, and 0=order parameter. 

On the other hand, the nonlinear interactions among atoms 
controlled by Feshbach resonance technique play an im- 
portant role for the ground-state properties. Fig. 3 plots a 
zero-temperature phase diagram for the atom-atom interac- 



tion strength v and the detuning A with a Rabi frequency 
in the framework of mean field. The Table lists the cor- 
responding ranges of the mean intracavity photon number 
(a'a\, the atomic population (J z ), and the "susceptibility" 
d (J z ) jdv for three different quantum phases. In the case of 
the repulsive interaction (v > 0), the critical point becomes 
A c = ±(u + v), which implies that an effective atom-field in- 
teraction is enhanced, while in the weak attractive interaction 
case (—it < v < 0), the effective interaction is suppressed. 
However, the basic features of the superradiant phases remain. 
In particular, in the case of v — —u, this system exhibits a 
novel second-order phase transition from the superradiant to 
the "Mott" phases (Red line) 025I1 . The relevant physics can be 
intuitively understood as following. In an optical cavity with 
co >> A, the cavity mode is only weakly or virtually excited, 



and the energy term oWa + (X/y/N) J x (aJ + a) is therefore 
nearly equal to — « J 2 /N. If v > — u, the ground-state proper- 
ties are governedby the energy — \u + v\ J x /N + AJ z + CtJ x . 
The effective potential in the Landau-Ginzburg theory is a 
double-well potential with the photon-assisted Josephson tun- 
neling, which means that this system is located at the su- 
perradiant phase. If —u — CI < v < —u, the energy 
— \u + v\ J z /N + AJ Z + ClJ x is dominant and the corre- 
sponding effective potential is a single-well potential with 
no internal Josephson tunneling, leading to the same atomic 
numbers for the two levels (m = 0), which may be referred 
to as the "Mott" phase J2^]. Also, when v is decreased, a 
second-order phase transition from the "Mott" to the "super- 
fluid" phases (Blue line) occurs at the critical point v c = 
—u — [CI 2 / 3 + A 2 / 3 ] 3 / 2 . In the so-called "superfluid" case, the 
effective potential is another double-well potential with the 
internal Josephson tunneling induced by the attractive inter- 
action i26ll . It should be pointed out that these three different 
phases can be distinguished experimentally by measuring the 
atomic population (J z ) and the "susceptibility" d(J z ) /dv. 
In the limit — > 0, this predicted second-order phase transi- 
tion from the superradiant to the "Mott" phases becomes a di- 
rect transition from the superradiant to the "superfluid" phases 
with the same order at the critical point v c = —u and A = 0. 

Although the second-order superradiant phase transition 
disappears in the strong attractive interaction (v < —u), 
another interesting phase transition (from the phase with 
nonzero macroscopic occupation of the level 1 to that of the 
level 2) in the "superfluid" regime emerges when the detuning 
A changes from negative to positive (i.e., from the red to the 
blue detunings). Fig. 4 shows the scaled atomic population 
m/N versus A for different Cls. We see that a novelfirst-order 
"superfluid" phase transition occurs at A = 0, and moreover 
this first-order phase transition exists until Cl c = \u + v\ (Red 
dashed line). For CI = Cl c , it becomes a second-order phase 
transition with the same critical point. For Cl > Cl c , no phase 
transition has been seen by varying A. 

We now estimate the energy scales for the parameters 
in Hamiltonian (2) to address the experimental feasibil- 
ity. Under the two-mode approximation, the wavefunc- 
tions of the macroscopic condensate states for the single 
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sion is anticipated 027I1 . 

In summery, we have established an extended Dicke model 
and designed a feasible experimental setup with controllable 
parameters. An exotic phase diagram has been obtained, 
which covers various phenomena from quantum optics to the 
BEC and reveals particularly several novel quantum phase 
transitions. 
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Research of China (No. 2006CB921800). 



FIG. 4: (Color online) The scaled ground-state atomic population 
m/N versus the detuning A for different Rabi frequencies Q when 

v < —it. 



magnetic trap may roughly be approximated by <pi(r) = 
n-z/^dxdydz)- 1 / 2 exp[-(x 2 /d 2 x + y 2 /d 2 y + z 2 /d 2 z )/2] with 



d x 



sJl/m R u) x , d y 



^1/mnujy and d z = ^Jl/m R uj z . 
Hence, the atom-atom interaction strength can be estimated 
by v = N(pi — pifi) / V^d x d v d z TTiji. For the typical val- 
ues (uj Xl ujy,Lu z ) =' 2tt x (290,43,277) Hz, p x = 4.2 nm, 
Pi, 2 = 9.7 nm, and mu = 1.45 x 10~ 25 kg, the energy scale 
of v is about -0.238 MHz with N = 5 x 10 4 , which en- 
sures that the error (the order of 1/ y/W) for determining the 
ground-state properties by means of the Holstein-Primakoff 
transformation is very low. The effective coupling st reng th 
A = 2Xy/N = 2.81 x 10 4 MHz for A = 2tt x 10 MHz Cj is 
indeed in the strong coupling regime. The e nerg y scale for u is 
about 0.315 MHz for w = 2.51 x 10 9 MHz Qj, which can be 
adjusted by controlling the frequency of photon. These energy 
scales for u and v imply that the intrinsic competition between 
the atom-atom and atom-field interaction should be taken into 
account seriously in the BEC coupled to the optical cavity. 
Also note that the aforementioned condition w < < u is well 
satisfied once w is tuned around Wx2 since u)\i ~ 6.8 x 10 3 
MHz << uj J23]. 

Finally, we elaborate briefly how to probe the predicted 
phase transitions experimentally. From the condition a — 
X(r] 2 - l)/2w(r/ 2 + 1) with A = 2.81 x 10 4 MHz and 
u> = 2.51 x 10 9 MHz, we can immediately evaluate the maxi- 
mum of the scaled mean intracavity photon number (a' a) /N 
and find it to be much less than the critical intracavity photon 
number n c = 7 2 /2<7q = 0.04. Therefore, one is able to per- 
form the transmission spectroscopy measurement with a weak 
probe laser to obtain the ground-state energy spectrum and 
atomic population since different quantum phases are, in gen- 
eral, characterized by their specific dispersion relations. The 
transmission (of this probe laser through the cavity) versus the 
detuning may be monitored and/or detected by counting pho- 
tons out of the cavity. Only when the probe laser frequency 
matches a system in resonance, the corresponding transmis- 
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